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The nature of the quantum valence transition is studied in the one-dimensional peri- 
odic Anderson model with Coulomb repulsion between f and conduction electrons by the 
density-matrix renormalization group method. It is found that the first-order valence tran- 
sition emerges with the quantum critical point and the crossover from the Kondo to the 
mixed-valence states is strongly stabilized by quantum fluctuation and electron correlation. 
It is found that the superconducting correlation is developed in the Kondo regime near the 
sharp valence increase. The origin of the superconductivity is ascribed to the development 
of the coherent motion of electrons with enhanced valence fluctuation, which results in the 
enhancement of the charge velocity, but not of the charge compressibility. Statements on the 
valence transition in connection with Ce metal and Ce compounds are given. 

KEYWORDS: CeCu2Ge2, valence transition, quantum critical point, valence fluctuation, su- 
perconductivity, DMRG, periodic Anderson model 

In CeCu2Ge2,^ the superconducting transition temperature Tsc has a maximum in the 
pressure regime far from the antiferromagnetic (AF) quantum critical point (QCP), where the 
coefficient A in the resistivity p = po + AT^ drops rapidly for n ~ 2. Since A is proportional 
to the effective mass^ (m*/m)'^ which is related to the f-electron number as m*/m = (1 — 
nf/2)/(l — nf),^ this implies that Tsc has a maximum around the sharp valence change of 
Ce. Similar behavior has been observed in the isostructural compound CeCu2Si2^ and in 
CeCu2(Sii_xGex)2-^ Recently, an NMR measurement has revealed that Tsc increases under 
hydrostatic pressure in Celrlus although AF spin fluctuation is suppressed,^ suggesting the 
different mechanism for the superconductivity from ordinary AF spin fluctuation. 

As known as the 7-a transition,^ Ce metal shows the first-order valence transition in 
its temperature-pressure phase diagram. The valence change of Ce is due to the 4f level lo- 
cated near the Fermi level, which can hybridize easily with the conduction band by applying 
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pressure. At the critical point, the valence susceptibility diverges as diverging density fluctu- 
ation in the liquid-gas transition. When the critical temperature is suppressed by controlling 
the chemical substitution and applying pressure, and enters the Fermi-degeneracy regime, 
diverging valence fluctuation is considered to be coupled with the Fermi-surface instability. 
This multiple instability in the quantum-degeneracy regime seems to be a key mechanism 
for understanding the instabilities observed in the above Ce compounds. However, theoretical 
understanding of the mechanism has not been fully achieved although a pioneering work on 
the valencc-fluctuation-mediated superconductivity^ and later discussions^ have been made. 

In this Letter, we clarify the nature of the first-order valence transition as well as the elec- 
tronic states near the quantum critical point. We show that the superconducting correlation 
becomes developed in the Kondo regime near the sharp valence crossover on the basis of the 
density- matrix renormalization group (DMRG)^*^ calculation for the one-dimensional periodic 
Anderson model 
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Here, the notation is standard and the number operator is defined by nf^ = al^cLia and 
nf = -I- for a = f and c. The last term is the Coulomb repulsion between f and 
conduction electrons, which is important in causing the valence transition.^' ^^'^^ 

The filling n is defined by n = (nf -|- nc)/2 with ria = J2iLi{''^i) / ^ for a = f and c with 
N being the number of lattice sites. We note that half filling is realized when n = 1 in this 
definition. In this Letter, we consider the case n = 7/8, t = 1, V = 0.1 and U = 100. In the 
DMRG calculation, we have kept the number of states up to 1500 and the system sizes up to 
N = 80 with the open-boundary condition. 

Figure 1(a) shows the f-electron-number nf vs the f-Ievel £f for U^c = 7.0, which has been 
obtained by a linear extrapolation to iV — oo using the data for N = 24, 32, 40, 48, 56, 64, 72 
and 80. We see that a jump in nf appears at £f = —5.6395, which indicates the first-order 
quantum phase transition, since a jump in nf results in the level crossing of the ground states 
by the relation nf = d{H)/d£i. The first-order valence transition (FOVT) between the Kondo 
state with nf 2± 1 and the mixed valence (MV) state with nf ~ 1/4 (= 2(1 — n) for 1/2 
< n < 1) is caused by Ufc, since a large Ufc forces the electrons to pour into either the f level 
or the conduction band. 

When Ufc decreases, the magnitude of the jump in nf decreases. We see that the jump 
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Fig. 1. Tif vs £f (open square) and Kp vs £f (solid diamond) for t = 1, V = 0.1 and U = 100 at 
n = 7/8 for (a) Ufc = 7.0, (b) 5.0 and (c) 0.0. 



is invisible at Ufc = 5.0 and the sharp crossover from the Kondo to MV states appears in 
Fig. 1(b). These results imply that the QCP at which both the jump disappears and the 
valence susceptibility Xf = —d'^{H)/def = —dn{/dei diverges, exists between U{c = 5.0 and 
7.0. As Ufc is set to be smaller than 5.0, the change in rif becomes smoother and finally a 
gradual crossover appears at Ufc = as seen in Fig. 1(c). 

To determine the location of the QCP, we have calculated the jump of the f-electron 



number, Arif 



ni 



n 



MV 
f ' 



at the FOVT point as a function of £f for several Ufc values. Here, 



the superscripts K and MV specify the Kondo and mixed-valence states, respectively. Then, 
the QCP is estimated at {e'^'^^ ,U^'~^^) ~ (—4.5,5.9) at which the jump disappears. In the 
slave-boson mean-field (MF) theory^"^ the QCP is identified at (e^^^, C/^^^p) = (0.19,0.98). 
The critical nature of the valence transition will be discussed in detail in a separate paper. 

We summarize the ground-state phase diagram in the C/fc-ef plane in Fig. 2. The solid line 
with solid diamonds represents the FOVT line determined by the DMRG, and the dashed 
line with open squares represents the point at which Xf has a maximum, which indicates the 
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crossover with enhanced valence fluctuation. Here, we note that the slope of the FOVT line 
is expressed as the ratio of the jump in rif and {n\n'i) at the transition by the 'Claudius- 
Clapcyron relation'/^ 6Ui^/5ei = -{nf - nf^)/{C^ - C}^'^), where Cfc is defined as Qc - 
Foi' ^fc > U^'^^ , 5Uic/Sei is shown to be in the MF theory/'' which 

coincides with the equation of the energy valance of the one-body f-electron energy at the 
FOVT; Ef + U{cn^ ~ ^, with fi being the chemical potential.^ However, it should be stressed 
that the correct slope for S> t^f^^^ is given as 6Uic/6e{ = — 1 by inserting the atomic- limit 
values of = 1, nf^^ = 1/4, = 3/4 and Cf^^ = into the thermodynamic relation. This 
discrepancy is due to the insufficiency of the description of the MV state by the MF theory. 

Here, we give the following rigorous remarks on the valence transition: For general filling 
ofl/2<n<l, the maximum value of the jump of nf, which is realized for large values of |ef| 
and C/fc, is given by Anf jmax = 2n — 1. This is in sharp contrast to the MF result, that always 
predicts the 4f^(Ce^+) to 4f''(Ce^+) transition; i.e., Anfjmax = 1- The slope of the FOVT line 
for large values of |ef| and Ufc is expressed by SUfc/Ssf = —1, irrespective of total filling nM 
It is noted that these results hold even in the square- and cubic-lattice systems. 

We also point out that the reason why the FOVT is difficult to realize in the Ce compounds 
can be understood by the smallncss of Ufc- In Cc metal, C/fc is large because of the on-site 4f-5d 
interaction, whereas in Ce compounds life is considered to be smaller even if the coordination 
number between Ce and surrounding atoms is taken into account, since U{c corresponds to 
the intcrsitc 4f-(d, or p) interaction. 

Wc have confirmed that at both sides of the FOVT, the peak structure appears at q = 7r/4 
in the Fourier component of {S^'^) for a=f and c when we apply a small magnetic field as a 
perturbation at the edge sites. The above period of the Friedel oscillations corresponds to 
2A;^^^ = vr/4, since kp"^^*^ = (1 -|- 3/4)7r/2 = Tvr/S. The large Fermi surfaces detected in the 
Kondo and MV states are consistent with the existence of the QCP, since by dctouring around 
the QCP, the Kondo and MV states can be adiabatically connected, with Luttinger's sum rule 
satisfied. 

In Fig. 2, the solid line with solid triangles represents the FOVT line determined by the 
slave boson MF theory. It is evident that the QCP (solid circle) identified by the DMRG is 
shifted from the QCP (open diamond) by the MF to the large-|ef | and large-C/fc direction. We 
also find no evidence of the phase separation at the FOVT, such as negative compressibility, 
and an inhomogeneous charge distribution in the DMRG results for the ground state. This 
result is due to the fact that the order parameter nf is not a conserving quantity, which is 
in sharp contrast to the MF theory where the FOVT is accompanied by the phase separa- 
tion.^^ These results indicate that the phase separation is destabilized by quantum fluctuation 
and electron correlation, and as a result the crossover regime from the Kondo to MV states is 
strongly stabilized. In connection to the experiments of the Ce compounds such as CeCu2Ge2, 
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Fig. 2. Ground-state phase diagram in C/fc-ef plane for t = 1, ^ = 0.1, U — 100 at n = 7/8. The first- 
order valence transition line (solid line with solid diamonds) with the quantum critical point (open 
circle) is determined by the DMRG. The dashed line with open squares represents the crossover 
determined by the point where Xf has its maximum value. The first-order valence transition line 
(solid line with solid triangles) with the QCP (open diamond) is determined by the slave-boson 
mean-field theory. 

this sharp crossover regime with strong fluctuation is most interesting, since the superconduc- 
tivity is observed in this regime. 

To clarify this point, we have calculated the superconducting correlation functions 
{P^±{i)P^J{i + x)) with P^±(i) = {ai^bi+rni ± ai|6i+m|)/\/2 for the singlet (-) and triplet 
(+) pairing. Here, a, b, c and d denote the f and/or conduction electron operators and we have 
calculated the nearest-neighbor (NN) pairing (m = 1) and also the onsite pairing (m = 0). 
To avoid the effect of the open boundary, we here fix the position of the operator at the 
central site, i.e., i = N/2 for the onsite pairing and i = N/2 — 1 for the NN pairing. Fig- 
ure 3 shows the superconducting correlations for Sf = —3.9 at Ufc = 5.0, which is close to the 
sharp valence change (see Fig. 1(b)). A remarkable point here is that the power of the decay 
of the correlations is smaller than 2 as clearly seen, for example, in the data of |(Pgl-Pfi^)| 
and K-Po^^ol^)! with dashed lines indicating x"^. We have confirmed that all correlations 
of \{P^_{{)P^_\i + x))\ show nearly the same power which is smaller than 2. On the other 
hand, for Sf = —4.2, which is in the Kondo regime, we find that the power of the decay of the 
pairing correlation is larger than 2. It is also found that the pairing correlation decays rapidly 
as ~ for Ef = —3.71, where the sharp valence crossover occurs. We confirmed that these 
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Fig. 3. Singlet pairing correlation functions \{Pm-i''')Pm- (* + ^))\ ^'^^^ Pm-^^) ^ {^i^hi+mi — 
arih+m^)/y/2 (see text) for ef = -3.9 and C/fc = 5.0 at t 1, F 0.1, U = 100 and n = 7/8: 
l^pccpcct^l (black open circle), |(Pf_Pf_^)| (black solid circle), \{P^IP1'l^)\ (pink diamond), 
|(pff pec t^)| (green triangle), \{P^lP§l'^)\ (orange triangle) and \ {Pi^PqJ)\ (blue square). Here, 
i = N /2 is set for m — Q and i — N/2 — 1 for m = 1. The dashed line ~ x"^ is drawn for 
comparison. 



results are not changed by the calculations at several system sizes. 

For U = Ufc = 0, the f level and conduction band are hybridized and /x is located at the 
lower hybridized band with kp^^*^. We confirmed that the charge and spin gaps are closed 
with kp"^^^ in the bulk limit in both the Kondo, MV and crossover regimes. Although the 
justification requires further careful studies, we assume the system being described by the 
single-band Tomonaga-Luttinger (TL) liquid in this Letter. For the single-component TL 
liquid without the spin and charge gaps, the singlet pairing correlation decays as x~^~^^^p }^ 
Since the charge and spin correlations decay as x~^~^f,^^ the above result implies that Kp > 1, 
namely, that the superconducting correlation becomes dominant. 

In order to estimate the TL parameter Kp systematically, we made a least-squares fit for 
the charge, spin and superconducting correlations defined by (OjOj) for Oi = Oi — (Oi) and 
Oi = nj, S^^ and P^_{i), respectively, since the conduction electrons have a large amplitude 
in the pairing correlation as seen in Fig. 3, which seems to give a reliable estimate. By assum- 
ing the asymptotic forms of the correlation functions of the single-component TL liquid, we 
obtained three -fCp's from the singlet-superconducting (SS), the charge-density- wave (CDW) 
and the spin-density-wave (SDW) correlation functions. Then, we take the middle point be- 
tween their maximum and minimum values as a mean value, with the error defined by the 
difference between the middle point and the maximum (minimum) value. 

The Kp values evaluated in this manner are shown in Fig. 1 as the solid diamonds. Here, 
for ef> — 3.75 at C/fc = 5.0 and for e{> — 5.6 at Ufc = 7.0, Kp is estimated from 4/cfCDW 
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and SS correlations. For Ufc = 5.0, Kp appears to exceed 1 at around Ef = —3.9 in Fig. 1(b), 
suggesting that the superconducting correlation is dominant, which is consistent with the 
result presented in Fig. 3. For Ufc = 7.0, Kp jumps at the FOVT as in Fig. 1(a). For Ufc = 0, 
Kp estimated from the SS correlation remains smaller than 1, indicating that the pairing 
correlation does not develop as in Fig. 1(c). These results show that the superconducting 
correlation is enhanced in the Kondo regime near the sharp valence crossover, which is close 
to the QCP of the valence transition. 

To get further insight into the nonmonotonic behavior Kp around Sf = —3.9 at Ufc = 5.0, 
we calculated the charge compressibility k = d{2n)/dfj,, since the TL parameter is expressed by 
the charge velocity Vc and k as Kp = ttVch/2. In Fig. 4, we plot the compressibility extrapolated 
to the N^oo limit in the form of k = 2/(A^Ac) with the charge gap Ac = [E{2nN + 2, 0) + 
E{2nN — 2, 0) — 2£'(2nA^, 0)]/2, where E{Nc,S) is the ground-state energy for the electron 
number A'c and total spin S. As ef increases, k, increases, since total charge fluctuation becomes 
large. We see that k does not show any diverging enhancement at the sharp valence crossover, 
Ef ~ — 3.71, in contrast to the enhanced Xf • This is due to the cancellation of the leading terms 
between (nffij) for a=f,c and (n[n^) in {{nl + n'i){n^j+nj)), since the diagonal and off-diagonal 
charge correlations have the opposite signs for Ufc > 0. 

A remarkable result here is that the compressibility stays constant around Ef = —3.9, 
where Kp has a maximum. This indicates that the enhancement of Kp is caused by the 
enhancement of the charge velocity, but not by the compressibility. This is quite different from 
the ordinary models, which exhibits Kp > 1.^^ In order to get further insight into this finding, 
we calculated the charge velocity defined by v* = limjv^oo ^c/^^/ with Aq = tt/{N + 1) 
instead of Vc,^'^ since the numerical accuracy is not enough for the second- lowest eigenvalue in 
the Lanczos diagonalization. We see that the maximum of v* appears around Ef = —3.9 and 
this suggests that the enhancement of is the origin of the enhancement of Kp. In Fig. 4, we 
show Kp = 'KKV*/2 for ef< — 3.85 by open diamonds and note that K* coincides with the Kp 
obtained by the least-squares fit of the correlation functions within the error bars. Since K* 
is obtained by v* and n after the extrapolation to AT — > oo, this suggests that Kp is evaluated 
properly within the finite sizes calculated in this study. A more accurate determination of Kp 
in the bulk limit is left for future study. 

The increase in v* toward ef~ — 3.9 in the Kondo regime is considered to be due to 
the decrease in the density of states at the Fermi level, i.e., the broadening of the Kondo 
peak. As ef further increases, is suppressed, since t/fc interrupts the coherent motion of 
the electrons in the MV state. This is the reason why has a maximum near the sharp 
valence increase. The enhanced implies the enlargement of the effective band width and 
thus the coherent motion of electrons with enhanced valence fluctuation is considered to be 
the origin of the development of the pairing correlation. The enhancement of Kp around 
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Fig. 4. Compressibility (solid triangle) obtained by k = liniAr^oo 2/(A^Ac) and charge velocity (solid 
circle) obtained by v* = liniAr^oo Ac/Aq for t ^ I, V = 0.1, U = 100 and f/fc = 5.0 at n = 7/8. 
The Kp (shaded diamond) obtained by the least-squares fit of correlation functions, K* = Trv*K/2 
(open diamond), and nf (open square) are shown. 

Ef = —3.9 also agrees with the following naive expectation: In the deep Kondo regime, the 
SDW correlation is dominant by the RKKY interaction and at the sharp valence crossover 
point the CDW correlation becomes dominant with enhanced Xf • Between the two regimes, the 
superconducting correlation can be dominant, since the regime midway between the Kondo 
state and the sharp valence crossover is unfavorable for both the SDW and CDW correlations. 

To summarize, we have clarified the novel nature of the electronic states near the QCP of 
the valence transition. Since valence fluctuations are basically ascribed to atomic origin, we 
believe that the basic nature of the valence transition and effect of its fluctuation revealed in 
this study is universal in two- and three-dimensional systems as well. 

One of the authors (S. W.) thanks H. Harima and S. Kawasaki for fruitful discussions. 
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